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Abstract 

The parity and time-reversal invariant effective lagrangian for a heavy fermion interacting 
with an abelian gauge field, i.e., NRQED, is constructed through order 1/M 4 . The imple- 
mentation of Lorentz invariance in the effective theory becomes nontrivial at this order, 
and a complete solution for Wilson coefficient constraints is obtained. Matching condi- 
tions in the one-fermion sector are presented in terms of form factors and two-photon 
matrix elements of the nucleon. The extension of NRQED to describe interactions of 
the heavy fermion with a light fermion is introduced. Sample applications are discussed; 
these include the computation of nuclear structure effects in atomic bound states, the 
J> model-independent analysis of radiative corrections to low-energy lepton-nucleon scat- 

tering, and the study of static electromagnetic properties of nucleons. 



1 Introduction 



Nonrelativistic QED (NRQED) is an effective field theory pQ describing the interactions of 
nonrelativistic fermions with electromagnetic fields. NRQED interactions at order 1/M 4 have 
become relevant for describing radiative corrections to proton structure contributions in hy- 
drogenic bound state spectroscopy [21 E] • The NRQED lagrangian, properly constrained by 
Lorentz invariance, trivializes the derivation of low-energy theorems of Compton scatter- 
ing (3] and automatically incorporates the intricate singularity structure of scattering am- 
plitudes [6]. It can be used to rigorously compute radiative corrections to low-energy 
lepton-nucleon scattering, and it also provides a model-independent framework within which 
to analyze static properties of nucleons, such as polarizabilities and generalized electromagnetic 
moments [7] . In this paper we derive a complete basis of operators and coefficient constraints 
through order 1/M 4 for the effective theory of nonrelativistic nucleons and leptons interacting 
with photonsj^] 

An important formal issue first arises at order 1/M 4 . As recently discussed in jH], a 
"reparameterization invariance" ansatz for enforcing relativistic invariance breaks down at 
this order. We derive the correct implementation of Lorentz invariance constraints and the 
resulting Wilson coefficient relations (i.e., nonrenormalization theorems) through order 1/M 4 . 
For applications involving NRQED of the proton or neutron we perform leading order (in a) 
matching computations to relate the remaining undetermined coefficients to observables of 
lepton-nucleon scattering. 

The remainder of the paper is structured as follows. In Section [2] we construct the NRQED 
lagrangian in the one-fermion sector through order 1/M 4 . This lagrangian describes the in- 
teraction of the proton with arbitrary background electromagnetic fields. Section [3] enforces 
constraints on the Wilson coefficients deriving from relativistic invariance, first by a variational 
calculation and then by an equivalent invariant operator construction. In Section [4] we perform 
matching calculations in the one-fermion sector, relating the undetermined Wilson coefficients 
to form factors and observables of (virtual) Compton scattering. In Section [5] we complete the 
basis of operators in the zero-fermion and two-fermion sectors, required to describe a proton 
interacting with a nonrelativistic lepton and dynamical photon. For applications to electron 
scattering (m e , E <C M) we also consider the extension to the case of a relativistic lepton. 
Section [6] provides a brief discussion of applications to the computation of nuclear structure 
effects in atomic bound states; to the model-independent analysis of radiative corrections to 
low-energy lepton-nucleon scattering; and to the study of static electromagnetic properties of 
nucleons. Section [7] provides a concluding discussion. 

2 Lagrangian 

Consider the lagrangian for a heavy fermion coupled to an abelian gauge field. We enforce 
hermiticity and invariance under parity, time-reversal and rotational symmetries. We also 
perform field redefinitions to eliminate time derivatives acting on the fermion field (apart 

1 For definiteness we will often refer to the heavy fermion ip as the "nucleon" , and to a second fermion \ 
in Section [5] as the "lepton" . 
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from the leading term); we refer to this choice as the "canonical form" of the heavy particle 
lagrangian. We thus find in the one-fermion sector through 1/M 4 , 

D 2 a B [d-E] a ■ (D x E — E x D) 



„ , t f.„ D 2 D 4 a B [d-E] 

C = * \ lDt + C2 2M + C4 8i^ + CF9 W + ^ W + lCs9 



8M 2 



{D 2 ,cr-B} DV ■ BD i cr DB D + D Bcr D 

+ ° Wl9 ^NP CW29 4M3 + Cp ' p9 8M^ 



{D\[d x B} 1 } 2 B 2 -E 2 2 E 2 



+ ""9 m + ca,9 - c A29 16M3 

\D 2 , D ■ E + E ■ D] {D 2 ,[d-E]\ \d 2 d ■ E] 

+ c xi9~ — ' vta ~ + C X29- -T71 — - + cxzg- 



M 4 M 4 A ™ M 4 

2 {D\ [E x B]*} DV ■ (D x E — E x D)D l e ijk a i D j [d ■ E]D k 
+ ^ Xi g w , + tc X5 g w , + ic xe g w , 

9 a-B[d-E] JE dcr B] JB dcr ■ E] 9 \E l a ■ dB 1 ] 

+ CX79 773 + C X89 77Z + °X99 77Z ~ + c xw9 



2 [BV • dE*] 2 a-Ex [d t E — d x B] 



+ c xll 9 2 1 M4 1 + c xl2 g 2 ^ i + G(l/M 5 )|^ • (1) 

We have defined D t = d/dt + igZA , D l = d/dx l — igZA 1 , where — gZ = — e, +e or 
for an electron, proton or neutron, respectively. We use the summation convention X % Y % = 
Yf i=1 X l Y\ and define [X, Y] = XY — YX, {X, Y} = XY + YX to denote commutators and 
anticommutators as usual. Square brackets around quantities imply that derivatives act only 
within the bracket. Electric and magnetic fields are defined as usual by E = — [d t A] — [dA°] 
and B = [d x A]. By the definition of E and B, [d ■ B\ = and [d t B + d x E\ = 0. 

The most general term in is obtained by constructing all possible rotationally invariant, 
hermitian combinations of iD t , D l , E*, iB l , ia l , with parity requiring an even number of factors 
of D l and E l . The operators through 1/M 3 were previously introduced in [TJ [§1 HD]- Terms 
at 1/M 4 with two field strength factors E % or B % are straightforward to tabulate; note that 
we have used [dtB] = — [d x E] and the assumption of canonical form to eliminate time 
derivatives of the magnetic field. Remaining terms at 1/M 4 involve one factor of electric field 
E % and three spatial derivatives D l . Spin-independent terms are straightforward to tabulate; 
the basis of operators parameterized by Cxi, Cx2, Cx3 differs from other possible choices by 
terms involving commutators [D\ i.e., terms with two field strengths. For spin-dependent 
terms we use [d x E] = — [dtB] and the assumption of canonical form to eliminate occurrences 
of [d x E]. The three- vector identity, 

D\E x a-y + (cr x D) j E i + a\D x E) j = D ■ E x a5 ij , (2) 

applied to remaining terms of the form ip^D l (. . . )D^ip, leaves the basis of operators parame- 
terized by Cx5, cxe- 
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3 Relativistic invariance 



The lagrangian ([I]) is invariant, by construction, under rotations and spacetime translations. 
The remaining constraints of relativity are enforced by demanding invariance under boosts. 



Here we derive these additional constraints, first by a variational calculation in Section 3.1 



and then by an equivalent invariant operator construction in Section 3.2 



3.1 Variational method 

As detailed in [8], under infinitesimal boosts, with infinitesimal boost parameter rj = —q/M, 
we may choose the heavy fermion to transform as 



2M 2 4M 4 

ic D g 



4M 2 



1 + 



D 2 

AM 2 



cs9 



+ ^-E+^q-<TxE + 0(g/M\l/MS) + ...j^, (3) 
while derivatives and gauge fields transform as Lorentz vectors: 

B^B-j^qxE, E^E + j^qxB, D^D+j^qD t , D t -> D t + j^q ■ D . 

(4) 

Field strength-dependent terms in pi) have been chosen to maintain canonical form. Since 
we are interested in the canonical lagrangian through order 1/M 4 , we need not specify the 
explicit form of the order 1/M 4 field strength-dependent terms, denoted by 0(g/M A ). A 
straightforward computation yields 



(5) 



where 

6fr = ^ [(1 - c 2 )iq ■ D] ip 
5C 2 = V t 

' 9 



--(l-c 2 ){q-D,D t } + ^(Z-2c F + c s )<T xq-E 



TP, 



sc 3 = ^ 



q-[8xB] (c F -c D + 2c M ) + -{q -D, D 2 } (c 2 - c 4 ) 



+ -f{q- D,a -B} (c 2 Z + 2c F - c s - 2c wi + 2c W2 ) 
o 

+ ^{<r • D, q ■ B} {-c 2 Z + c F - c p/p ) + ^q cr(D B + B D) (-c F + c s 



VP i 



'(6) 
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From 6C lt SC 2 and 5C 3 , we find PH E0 

c 2 = 1 , c s = 2c F - Z , c 4 = 1 , 2c M = Co - c F , c W2 = c wl - Z , c p / p = c F - Z . 

(7) 

Employing the above relations, the variation 5C4 takes the form 

5£ 4 = ¥ j [D 2 ,q ■ E] (^- Cf + c d - 32c x ^J + ^{q D,[d E}} (~ + c F - 16c X2 

a 2 I ' Z 2 
+ y —q ■ E x B [ — + 2c F {Z - c F ) - 2Zc D + c A2 + 16c X4 

+ ^[q-(T x dd ■ E] (-Z + c F - ^c D + c wl + 8c X6 



+ |l7 (q\E x <r) J ' + x <r)V + a x q ■ E5 tj ) D j - 2c F + 16c X5 ^ 



if, 



where we have suppressed terms that are removed by field strength-dependent modifications 
of the boost generator, denoted by 0(g/M 4 ) in ([3]). We readily find, 



32cxi 


5Z 

~ T ' 


- c F + c D , 


32c X2 


Z 
~2 


+ 2c F} 


32c X 4 


= -Z 2 


- Ac F (Z - c F ) +AZc D 


32cx5 


= -Z 


+ 4c F , 


32c X6 


= 4(Z 


- c F ) + c D - Ac W i , 



(9) 

while coefficients cxz an d c X7...xn are n °t constrained by Lorentz invariance. We thus find 
that seven new quantities are required at order 1/M 4 to describe the proton's response to 
arbitrary background electromagnetic fields. 

3.2 Invariant operators 

An alternate method for enforcing Lorentz invariance is to construct the lagrangian from 
explicitly invariant operators. We summarize here the main points; for details see [H]. 

The basic building block in the construction is the field = T(v,iD)ip v , where ip v is a 
Dirac spinor field with ^tp v = ip v . The matrix- valued operator F(v,iD) is determined such 
that under an infinitesimal boost A, where K^ v v v = v M + q^/M, the field \E% has a simple 
transformation law: * v ->■ e iq ' x ^ v . Noting that e'^iiD 1 " + Mv^ + q»)e iq ' x = iD» + M/, 
we may thus build invariant bilinears from contractions of polynomials of 7^ and V M = + 
iD^/M, between % and \EV 



As noted in [3J, we find the opposite sign in the relation for cm in nfl compared to [10] . 
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The function T(v,iD) is a solution to the invariance equation, 

T(v + q/M, iD - q)A~ 1 W(A } iD + Mv) = T(v, iD), (10) 

where W(A,p) is an element of the little group for timelike invariant vector v^, following 
from the theory of induced representations of the Lorentz group. Up to the relevant order for 
determining the 1/M 4 lagrangian we have [8] 

r = + J_ i--UD ± ) 2 - -ip ± iv ■ d] + J- [-UDAHv ■ D 

2M M 2 \ 8 V ; 2 r I M 3 l4 V ±; 



~(iD ± ) 2 + (ivD) 2 



+ g -fF, v v»Dl + ^a^F^ip j.} + • • • , (11) 



where we have defined Z)^ = -D M — i^u • -D, and for abelian gauge fields F M „ = d^A u — dyA^. 
Note that the last two terms of (11) are absent in the ansatz for reparameterization invariance 
given in [11], leading to incorrect Lorentz invariance constraints at 1/M 4 and beyond. This 
subtlety is explained in [8]. 

A complete basis of invariant bilinears required through order 1/M 4 is 

o^Fyv , . {V„[MV U ,F^]} [MV a ,[MV a ,^F F ]] 



C = ty v < M(Y — 1) — a F g + iaog — awi9~ 



AM ' ^ 16M 2 VVL " 16M 3 

12 



+ + a ^ 2 Va 7 Q M^ h + axsBx3 + I>^" ( 12 ) 

The bilinears Bxi for i = 3, 7 ... 12 are chosen to reduce to the respective operators multiplying 
Cxi in upon setting = (1,0,0,0) and neglecting 1/M suppressed corrections. Since we 
are concerned only with the lagrangian through order 1/M 4 we do not specify an explicit 
choice for these Bxi- A computation shows that the field redefinition to recover canonical 
form is 

1 s^^fi iv-D\ gZ gZ igZ n u ™„ 



8M 3 1 4M 3 x -^-l""' J 8M s 

+ W a - [ ^'^ ] }^- (13) 

Upon setting = (1,0,0,0), the resulting lagrangian, expressed in terms of ip' v , is identical 
to the previous result with constraints ^ and 



4 Matching: one-fermion sector 

In contrast to NRQED for the electron or other fundamental fermions, the matching for a 
composite particle such as the proton cannot be performed perturbatively. We instead must 
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appeal to nonperturbative (e.g. lattice) methods, or to experimental measurements. This 
section relates the matching conditions in the one-fermion sector to standard form factors and 
two-photon matrix elements of the nucleon. 



4.1 One-photon matching 





Figure 1: Tree level matching of the one-photon amplitude in the full theory and NRQED. 
The black dot in the diagram on the right-hand side represents single-photon NRQED vertices. 

Consider first the operators contributing to the one-photon matrix element of the nucleon. 
The matching is performed in terms of standard invariant form factors, 



i a 



(N(jpf)\Jjr\N(p)) =v{p')VM)<P) , r„(g) = 7m W) + ^W)?* 



(14) 



where q = p' — p and N denotes a proton or neutron; we suppress the superscript N in the 
following. Equating the effective theory with the full theory^ we find (cf. Fig. [T]) 

c F = Fi + F 2 = Z + a N + O{o) , 

c D = F 1 + 2F 2 + 8F[ = Z + ^M\r N E f + 0(a) , 



cwi = F 1 + -F 2 + AF[ + AF' 2 
cxs = \f[ + \f> 2 + \ F '( , 



(15) 



where Z denotes the electric charge, a at is the anomalous magnetic moment of the nucleon, 
and is the nucleon charge radius. We have introduced dimensionless barred quantities to 
denote derivatives with respect to q 2 /M 2 at q 2 = 0: F\ = Fi(0) = Z, F 2 = F 2 (0) = a N , 
F[ = M 2 F[(0), etc. The new quantity F" appears at 1/M 4 . Expressions for other Wilson 
coefficients through 1/M 3 in terms of form factors can be found using 0. At 1/M 4 , we also 
find 

cxi = —F x + — F 2 + -Ft , 
X1 128 32 2 4 1 ' 



The nonrelativistic normalization of states in NRQED is obtained using u(p)u(p) = M/E p in (14) 
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CX2 



CX5 



lp x + -F 2 , 
64 1 16 ' 



3 

32 



—A 
32 



8 4 1 



2 2 



(16) 



and it is readily verified that these expressions satisfy the constraints ([9]). In the presence 
of radiative corrections, the form factors on the right hand sides of (15) and (16) should be 



interpreted in an appropriate infrared regularization scheme; alternatively, the matching may 
be performed with infrared finite observables. The corresponding infrared subtractions and 
ultraviolet renormalizations must be performed to obtain the Wilson coefficients including 
radiative corrections^ 



4.2 Two-photon matching 





p + v 





V + 




Figure 2: Tree- level matching of the Compton scattering amplitude in the full theory and 
NRQED. The black vertices in the diagrams on the right-hand side represent NRQED vertices. 

The Compton scattering process, 7*(<?)A r (p) — > r ){q l )N{p') ) with one virtual and one real 
photon is sufficient to determine the remaining coefficients in the 1/M 4 NRQED lagrangian. 
Consider the low-energy expansion of the virtual Compton scattering amplitude as depicted 
in Fig. [2j Let us define a conventional separation 

jvi _ JVl Bmn -|- y^l non _Born ' \ l ' ) 

by declaring that "Born" terms are defined by the "Sticking In Form Factors" (SIFF) pre- 



scription using the form factors of (14). Explicitly 



1 



-M 



H-g') Up) 



18) 



The expressions on the right hand side of (15) and (16) correspond to those referred to as c^ ED in [9]. 



The renormalization procedure in dimensional regularization is described in [II 
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This convention ensures that q^M-B OVVL = ol-M-Bom = 0' so that the same Ward identities, 
^ M Zn-Bom = ^-^non-Bom = °; ma y be applied to constrain -A/Cn-Bom- We adopt a conve- 
nient basis for A^ n _ Born as given by Drechsel et al (Eq. (A10) of [T2]). 
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-^non-Born = Yl 9 ' 9 '' 9 ' P ^ " ( 19 ) 



1=1 



Subtraction of the Born terms ensures that the residual contributions fi(q 2 ,q ■ q',q ■ p) are 
free of kinematic singularities, and may thus be Taylor expanded at small photon energy. 
Employing q ■ p ~ oj, q 2 ~ u 2 , q ■ q' ~ w 2 , to the relevant order, we require 

U{ q \ q-q',q-p) = ho + 0(u 2 ) , « = i,2, 5, 6, 11, 12 , 
/i(g 2 ,g-g / ,g-p)=Aig-^ + 0(a; 3 ), z = 4,10, (20) 



where we employ the notation of 

The matching can be performed in arbitrary reference frame, with the results]^] 

4M 3 

c A1 = F 2 + 4/ 1)0 = Z 2 + p M + 0(a) , 

a 

CA2 = 4FxF 2 + 2F| + 16A^ + 32/ 2 , = 2c& + —M 2 {r N E ) 2 (a E + M ) + 0(a) 

S a 

cx7 = -^FiFl + 2^1-^2 ~~ 2^ 5>0 ~~ ^ f 11,0 — 2/12,0 , 

0x8 = ~§2 P " ~ l Pip2 ~ h p2 ~ ( Pi + l p2 ) p ' ~ \ Pip2 ~ 4/4,1 + 4/6,0 " 2/10,1 ' 



cx9 = + ^F X F 2 + l -F 2 2 + i(A + F 2 )(F 1 / + - 2/ 10> i , 

cxio = ^F 2 + ^F 1 F 2 + + 4/ 4>1 + 2/ 10>1 + 2/ u , , 

cxii = + ^AA - i(A + F 2 )(A + F2) + 2/ 10 ,i + 2/u,o , 

64 Id 2 



Lp 2 - ^F 1 F 2 - ^F 2 -\(F 1 + F 2 )F[ - + 4/ 4)1 + i/ 5>0 + 2/ 10 ,i + 2/ n ,o , 

(21) 



where the dimensionless form factor derivatives, F- n \ have been introduced after (15) and we 
have similarly defined dimensionless quantities f\ = M 3 fi, f 2 = M 5 f 2 , / 4 = M 4 / 4 , f 5 = M A f 5 , 
f e = M 5 / 6 , f w = M 3 / 10 , fix = M 4 f u , fi2 = M 5 f 12 . We have then used the expansion of the 
amplitudes fi given in (20), with dimensionless quantities = fi$+0{uj 2 ) (i = 1, 2, 5, 6, 11, 12) 

5 We have performed the computation both in the laboratory frame, p — 0, and in the center of mass frame, 
q + p = 0, by matching both the full theory and the effective theory to the twelve independent spin structures 
for the virtual Compton scattering process. 
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and £ = /; )ig - p/M 2 + CV 



4, 10). Note that six new phenomenological parameters are 



required at 1/M 4 to match cxi...x\2- An expression for cxa in terms of scattering observables 
can be found using (§. The expressions (plT) can be translated to a multitude of other 
conventions for the observables of Compton scattering]^] 



5 Pure photon and four-fermion operators 

So far our analysis has focused on the one-fermion sector. We have derived the form of the 
lagrangian appropriate, e.g., to a proton in a background electromagnetic field. Let us consider 
the complete QED theory including dynamical photon, as well as a lepton (electron or muon) 
field. The case of a nonrelativistic lepton is appropriate to bound state hydrogen studies, or 
very low-energy lepton-nucleon (e.g. muon-proton) scattering, where E <C M x , M. We first 
consider this case, constructing the operator basis, deriving coefficient relations and identifying 
redundant operators. We then turn to a brief discussion of the case of a relativistic lepton, 
appropriate to e.g. low-energy electron-proton scattering with me, E <C M. 



5.1 Pure photon operators 



The pure gauge sector for NRQED is the well known Euler-Heisenberg lagrangian. Enforcing 
parity and time reversal symmetry and neglecting total derivatives we find 

£ 7 - —F^F +c V2 — +c V4 — +c E1 g M4 +CE29 W + • • • • 

(22) 

The coefficients cy 2 and cva may be set to zero through field redefinitions on A M , as discussed 
in Section l5\3l below. 



5.2 Four-fermion operators 

Consider four-fermion operators relevant for processes in the one-nucleon, one-lepton sector. 
We enforce hermiticity and invariance under parity, time-reversal and rotational symmetries. 
We use the notation D for a covariant derivative acting to the left, X*D % = [d l X\ + igZXA 1 , 
and define D + = D + Having performed field redefinitions to eliminate 

operators with time derivatives acting on heavy fermions, the lagrangian in this sector, through 
1/M 4 , is 

= X j ^ X + jj^ty X ] X + X ] D\ X + ^&J> X j Dl X 



+ ^(D 2 + £v x ] x + x\n 2 + 



6 For example, Compton scattering of real photons determines all coefficients, apart from cx7 and Cxi2, in 
terms of the conventional electric and magnetic polarizabilities a e , (3m and Ragusa's [4 spin polarizabilities 
7l : M 3 a E /a_= -/y, - 4/ 2)0 , M 3 /3 M /a_ = f lj0 , M^/a = -8/4,1 - 4/ 10 ,i - 4/ Uj0j M^/a = 4/ 10ll , 
M i ^/a = 4/ 6 o + 2/n , M 4 7 4 /a = — 4/ 10 i — 2/n , where a is the fine structure constant. Compare with 
Eq. (28) of H3J'. 
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M 4 
d\ 
M 4 



M 4 ' 
M 4 



id 9 
M 4 ' 
id\ 
M 4 



X 



_o t^os do 



+ xV(£> 2 + S 2 ) X + ^Vty *V(LW + frfr) X • (23) 

Here x is the nonrelativistic lepton field with mass M x and for notational simplicity we write 
all operators in terms of the common mass scale Covariant derivatives appearing within 
a fermion bilinear in (23) are understood to act only on fields in that bilinear. The heavy field 
Tp transforms under boosts as in (J3|. Recalling that q in (|3| is related to the mass-independent 
infinitesimal boost parameter by r/ = —q/M, the transformation law for x is obtained by the 
replacement M — > rM and q — > rq, where we define r = M x /M. We thus find 



5C 



iffUq ■ D_ip x ] X 



M 4 

+ iftcr ■ q x D + 4> x j X 
+ ^Hx^° ■ Q x D+x 
+ iftUcr ■ D_ip x}cr " QX 

+ ipHcr ■ qD l _ip x}°' l X 
+ £>(1/M 5 ) . 



do 



2rd, - 2d 



d 2 
2r 



do d\ 
4 4r 

do di 
~ + ^-2rd n -2d 12 
4r 4 



d\ 



2rd 16 - d 19 



+ ipH x ] iq- D-x 

+ ipHq ■ D_a l ip x?a % x 
+ ip'cr ■ qip x i* 7 ' D-X 



2c?4 — 2rd§ 



d\ 
~2 

d\ 

2r 



— 2rdi4 — 2d\s 

- 2d u - 2rd 2 o 



di 
4r 

d\ 



2d 16 - rd 21 



—. 2d 17 - rd 2 i 

4r 



(24) 



This enforces the relations 

rd 4 + d 5 = ^-, d 5 = r 2 d 6 , 8r(d 8 + rd 9 ) = ~rd 2 + d t , 8r(rdu + d X2 ) = — d 2 + rd 1 , 



rd u + dis = % , dis = r 2 d 20 , 2rd 16 + d 19 = ^ 



r((ii 6 + dn) + ri 19 = , d\ 9 = r 2 d 



21 j 

(25) 



7 Note that the coefficients c?i.2 in (23) are related to those of Caswell and Lepage [I] by a factor M x /M. 
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implying a total of twelve independent four-fermion operators through 1/M 4 , including two 
at order 1/M 2 . By performing field redefinitions on the gauge field A^, some of these four- 
fermion operators are found to mix with one-heavy particle sector operators, as discussed in 
Section 15.31 below. 



The lagrangian (23), with constraints (25), applies to the case of distinct heavy parti- 
cles represented by ip, \i with arbitrary mass ratio M x /M. For certain applications, e.g. 
positronium or heavy quarkonium bound states, the fields if) and x can be taken to represent 
particle-antiparticle pairs with r = M x /M = 1. Charge conjugation symmetry is then im- 
plemented by enforcing invariance under if) ■<->■ x, thus reducing the basis of operators. This 
case has been investigated for QCD through 0(1/M 4 ) by Brambilla et al. [H]. We find that 



our basis of four-fermion operators (23) and constraints (25) are equivalent to those found in 
Ref. [Hj for this special casej^] 

5.3 Field redefinitions and redundant operators 

With a dynamical photon field, we may perform field redefinitions that maintain reality and 
gauge, parity, time reversal and rotational symmetries. In order to avoid upsetting the pre- 
viously determined coefficient relations, we must also maintain the transformation law for 
four- vector under Lorentz transformations, i.e., 

A°^A°-^-q-A, A^A-^-qA°. (26) 

Let us write 

A^ = ^ + A y A^ + A^ + A X A^ + ... . (27) 
For the pure gauge field terms the most general expression is 

A^ = 0,1-^5- + + 0(l/M G ) . (28) 

Terms involving the heavy fermion if) take the form 

a V>i — iTTo r a^ 2 h a^g 



g V1 M 2 WL M 3 ^ M 4 w * M 4 

v v a>«*{ve,F a p}* v 6 

+ a ^g j^i + 0(1/ M ) , (29) 



where we have employed the invariant operator formalism of Section 3.2 In particular, 



Tijj v with T from (11 ) and ifj v from (13), expressed in terms of the field ip' v = if) with canonical 
lagrangian ([I]). As an alternative to the invariant operator formalism employed in (29) we 
may expand A^A° and A^A in a series of rotationally invariant operators with arbitrary 



8 The difference between abeiian and nonabelian gauge fields is trivial for four-fermion operators through 
this order. 
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coefficients, and subsequently constrain these coefficients using (26). The result is equivalent 
to (29), with five free parameters through 0(1/M 4 ), 



9 
9 



CLipl 



-Q>i/)l- 



V 4 



2M 3 



+ a^p4 



M 3 



+ di>59 



ib^cr x Dib ib^cr ■ Bib m , , ^ 

- — w — + a ^ 9 ~^- + ° {1/M) 

ip^cr x Eib 



M 4 



+ C(l/M 5 ) . 



(30) 



The expansion of A x v4 M is obtained from (30) with the replacements ip — > Xi M — > M. 
Z — ^ Z x and a^i — > a X i. In terms of the field A' in (27), we find in the pure photon sector 

1 11 



x< 



5cv2 = — 2 fl 7i ) $ C V4: = ~2 a " /2 ~~ 4 a T 1 ^ ^ a "7 lCv2 ' 
while for the ib sector, 

5cd = —8Za^i + 8a^,i , Scwi = — 4cpa 7 i + 8a^,4 , 5ca2 = — 16Z 2 a 7 i + VoZa^\ , 
cdQ 7 i „ . r. csZa^i 



(31) 



-.2 



4 



5c 



X12 



Cp Z Clyi 

csZa^i 



2~~ 



— Za^ 4 , 5cx9 



h CpCt^ , OCxil 



Similar relations hold for the Wilson coefficients in the x lagrangian, defined as in (1), 
with ib — > x, Z — > Z x , M — > M x , Ci — > c\ . Finally, for the four-fermion operator coefficients, 



+ <^3 



CpCL-yl 



(32) 



5d 2 
9 2 ^ 


-Z x a^t - 


Za xl 


5^3 a^i cdci x i Za X 2 

g 2 ~ g r 2 1 g r 2 1 ~X°^ 1 r 4 ' 




Z x a^i 
Ar 


Za xl 

4 r 3 


5dj 


ZZ 


5d 8 

IF = 


— [a^i - 


- 4a^ 4 ) 


c s a x i 
8r 2 ' 


^2 = 4r 4 (°Xl 4 M r 4 , 


5d n 
9 2 


Cg Ct-ipl 

8r 2 


Z 

+ 8H ( 


a x i - 4a x4 ; 


g2 g2 2r 2r 3 



(33) 

The coefficient relations ([7]), ^ and (25) are preserved, since by construction the Lorentz 
transformation properties of A^ are unchanged and hence the boost transformation rules ^ 
and Q still apply. 

We may use (31) to eliminate vacuum polarization terms Cyi and cy^ in favor of compen- 
sating terms in (32). Similarly, (32), together with the analogous relations for c[ x \ and (33), 



can be used to eliminate 10 linear combinations of Wilson coefficients for two-fermion and 
four-fermion operators. Different applications may favor elimination of different operatorsr' 



We have not specified gauge fixing and source terms, which are also affected by field redefinitions. 
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5.4 Relativist ic lepton 

For applications such as lepton-nucleon scattering at energies mi,E <C M (e.g., low-energy 
electron-proton scattering), the relevant effective theory involves a heavy fermion (e.g., the 
proton) interacting with an electromagnetically charged relativistic fermion (e.g., the electron). 
Let us briefly discuss this case. Enforcing parity, time-reversal, gauge, Lorentz, as well as chiral 
symmetry at mi = 0, we find the leptonic interactions with the photon, 



Ct 



mi + gey m t ^ + gc\ ™i^ + gey — — 2 + 0(1/M 



(34) 



where we assume field redefinitions have been performed to remove power suppressed terms 
involving (ip — mi)l. 

Having performed field redefinitions to eliminate operators with time derivatives acting on 
fermion fields, the lagrangian for the nucleon-relativistic lepton sector through 0(1/ M 3 ) is 

Am = jj^V *7°* + Hid + m e ££ + ^iD^ Hi 

+ ^VW«7 • + ^e !J '^VV m e £a jk £ + ^e ijk ^a^ £^D k + £ 

+ ^V^V h°lsiDH + ^L^aHDii; £ 7 ° l5 £ + 0(1/M 4 ), (35) 

where £ is the relativistic lepton field with mass mi and a 1 ^ = ffyST 3 ']- The heavy field ip 
transforms under boosts as in (|3]), while £ transforms under finite dimensional representations 
of the Lorentz group in the usual way. Under Lorentz transformation, we thus find 

SCpt = -^H£q l l l £ (h + 26 4 ) - ^ W#rV (b 2 + 2b 9 ) + 0(1/M 4 ). (36) 

This enforces the relations ^ ^ 

h=~bi, b 9 = --b 2 , (37) 

leaving seven operators in this sector through order 1/M 3 , including two at order 1/M 2 . 

By performing field redefinitions on the gauge field A^, some of these four-fermion operators 
are found to mix with one-heavy particle operators. In addition to the contributions A 1 A tl 
and A^v4 M from (27) we may employ 

A e A» = gaJ-^ + 0(l/M*). (38) 

We thus find the modified couplings in £ e , 

8c { § = -Zia-fi + a El , (39) 
and for the four fermion operators in £^i, 

= —Zan — Zga^i , = —Zia^i — -cfclh , (40) 

g 2 g 2 2 



with relation (37) remaining intact. 
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6 Applications 



Applications of the NRQED lagrangian for the nucleoli include: the computation of pro- 
ton structure effects in atomic bound states; the model-independent analysis of radiative 
corrections to low-energy lepton-nucleon scattering; and the study of static electromagnetic 
properties of nucleons. Let us discuss sample applications in each of these areas. 



6.1 Proton structure in atomic bound states 

As a sample computation involving the 1/M 4 NRQED lagrangian J^j let us analyze the effects 
of nuclear structure on the two-photon exchange contribution to the 25 — 2P Lamb shift in hy- 
drogenic bound states; this contribution is the subject of intense scrutiny [16 ] \T7 \ \18 \ \19 \ 13] 120] 
due to the discrepant measurements in muonic and electronic hydrogen [2]. The dominant 
theoretical uncertainties in the muonic hydrogen Lamb shift arise from proton structure, repre- 
sented by contributions through 0{a) to cp in Q, arising from first order vertex corrections, 



and 0(a 2 ) contributions to d,2 in (23) arising from two-photon exchange [3J. We let if) of 
Section [2] denote the proton of mass M, and x °f Section [5] denote the electron (or muon) of 
mass M x = m e . Electric charge assignments are given after (JlJ. We focus here on the model- 
independent result for the leading two-photon exchange corrections in the limit m e /M <C 1. 
Structure dependent corrections to d<i lead to first order energy shifts in hydrogen p] 

tEM = (41) 



Tin 3 \ M 2 y 

where (n,£) are principal and orbital quantum numbers and m r = m e M/(m e + M) is the 
reduced mass. The matching condition for involves a weighted moment of the structure 
functions for forward Compton scattering on the proton. From [3J, the necessary hadronic 
input is contained in 

g = Jmg f dxV r^ r + + subtractions , (42) 

where the subtractions depend on infrared and ultraviolet regulators. The structure functions 
in this expression are evaluated at v 2 = —Ax 2 M 2 Q 2 . They are defined by the forward proton 
matrix element, 

W»(k,q,8) = i J d*se^<MT{J^(aO^(0)}|M>- ( 4 3) 



10 For bound state applications it may be computationally efficient to further integrate out offshell momen- 
tum modes and/or higher Fock states to arrive at an explicitly v/c expanded NRQED [HE], or fixed particle 
number quantum mechanics |15| . Our purpose here is to examine the impact of nucleon structure, as described 



by finite shifts in the Wilson c oeffi cients of the theory described by (|T|) , (p2| , ((23 1 
11 As mentioned in Section 



5.2 



we have used powers of 1/M (not 1/M X ) to define four-fermion Wilson 
coefficients in (23). We thus have dx^/M 2 — df^/(m e M) where d^ 1 " is the coefficient of Caswell and Lepage [J, 
also employed in [3]. 
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For the present application we require the spin-averaged component symmetric in /i, v, 



k ■ qq h 



W 2 (u,Q 2 ), 
(44) 



where Q 2 = —q 2 and v = 2k ■ q. Our normalizations are such that for a point particle, 
Wx = 2z/ 2 /(Q 4 - v 2 ) and W 2 = 8Q 2 /(Q 4 - 

Let us proceed to evaluate the finite shift 5d 2 due to proton structure, and hence 5E in 
(41), in the limit m e -C M. Details of this computation are presented elsewhere. The essential 
input from NRQED is the evaluation of W;(0,Q 2 ) through 0(Q 2 ), 

Wt{0, Q 2 ) = -2 + 2c 2 F + ^ (4 - 2c F c wl + 2c M + c Al ) + 0(Q 4 ) , 



M 2 W 2 (0,Q 2 ) 



8M 2 9 Q 2 

+ 32(C X 1 + C X 2 + CX3 



1 + c 2 F - c D + -c 2 D + 2c M - 2c F c W i 



-CA2 



0(Q 4 



(45) 



and the evaluation of Wi(u, 0) through 0{u 2 



Wx{v,$) = -2 + 



lim W 2 (z/,Q 2 ) = 



8M 4 



-c| + c F c s + 2c M - -C A 2 + 0{v 



(46) 



where the final relation follows from requiring that the Compton amplitude (44) is finite at 
Q 2 ->■ 0. Note in particular that W 2 (0,Q 2 ) relies on C(l/M 4 ) one-photon vertices involving 
Cxi, cx 2 , Cx3- This occurs because the 0{Q 2 ) terms in W 2 (0,Q 2 ) arise at third order in the 
Q 2 /M 2 expansion. Employing these results in the matching relation (42) yields the leading 
contribution to S'-state energies at m e /M — > 0, 



SE(nS) = log [M 3 (5a E - fa) /a - 3a p (l + a p ) + 2M 2 (r^) 2 ] +...}. (47) 

7T77. M 6 Ad 1 ' 

Our result for the coefficient of the logarithm containing a F and /3m was obtained in a temporal 
gauge analysis in [21] . We have here used NRQED to establish the complete logarithmically 
enhanced contribution. The expansion at small lepton mass is not possible for muonic hydro- 
gen owing to low hadronic scales scales m M ~ m n ~ — m p . Further analysis for this case 
will be presented elsewhere. 



6.2 Radiative corrections to lepton-nucleon scattering 

Extraction of key hadronic quantities from elastic lepton-nucleon scattering, such as the charge 
and magnetic radii of the proton, demands rigorous control over radiative corrections. Tradi- 
tional analyses of the two-photon exchange contribution to elastic electron-proton scattering 
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Figure 3: Diagrams contributing at leading order in 1/M to two-photon exchange in electron- 
proton scattering. 



resort to hadronic models involving phenomenological form factor insertions into point-particle 
Feynman diagrams [22| |23| EU [25] , or require modeling of soft functions to evaluate expres- 
sions obtained in a hard-scattering factorization framework [26J. At low energies, E <C M, 
a systematic and model- independent approach is provided by NRQED. This region has over- 
lap with current and planned electron-proton and muon-proton scattering measurements, and 
provides a rigorous test of hadronic models employed at higher energy. 

Consider the two-photon exchange corrections to low-energy electron-proton scattering in 
the limit rri£ = m e <C E <C M, where E is the energy of the incident electron in the rest 
frame of the initial-state proton. The appropriate effective theory consists of a relativistic 
lepton and heavy proton, as described in Section 5.4| The leading order NRQED diagrams 



are not sensitive to proton structure, yielding, in Feynman gauge, [electric charge assignments 
are given after Q] 

iM « e A u(p )7°7 M 7°u(p) 

d A L L„ 1 1/1 1 

x 



(2tt) 4 L 2 + iO (L - p) 2 - X 2 + iO (L - p') 2 - A 2 + iO V-^° + ^ + ^0 L° - E + iO 

(48) 

where we have neglected the electron mass and used that p° — p'°, k° — M, and k'° — M 
are sub leading in the 1/M expansion to simplify the integrand. A photon mass A is used 
to regulate infrared divergences. Evaluating the integral in the limit A E, we recover the 
expression for relativistic lepton scattering in an external Coulomb field [2T| [28J q 

Ana . . n . , f In Q ( A Q 2 C 

M = + + i ( - 2 log — + log l E 

+ 0[a 2 ,\/E,m e /E,E/M] \ . (49) 



12 There is an apparent typographical error in Eq. (2.6) of [28 where in the second line the integral / should 
have the opposite overall sign. This affects only the imaginary part of the 0(a) correction to the leading 
amplitude and is thus not relevant to first order radiative corrections in the cross section. 
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The power of the effective theory lies in the possibility to systematically compute corrections 
in powers of E/M and to relate observables such as scattering amplitudes and bound state 
energies, in a model-independent fashion. Sensitivity to the Wilson coefficient cf appears at 
order 1/M in the proton spin-dependent amplitude. At order 1/M 2 there is a dependence on 
Co, and on the constants 61 and b 2 from (35), which encode information on proton excitations 
and may be related to moments of the forward Compton amplitude. 

A similar analysis may be performed to compute radiative corrections to muon-proton 
scattering in the limit E <C ~ M, where now the effective theory consists of heavy muon 
and proton fields, as discussed in Section |5.2| Proton structure not captured by elastic form 
factors is first encountered at 0(1/M 2 ), encoded in coefficients di and d 2 in (23). NRQED 
may similarly be used to compute radiative backgrounds to searches for possible new low-mass 
and weakly coupled particles using low-energy electron-proton scattering 



6.3 Spin polarizabilities and static properties of nucleons 

Having constructed the NRQED lagrangian, it is straightforward to relate parameters ex- 
tracted from scattering measurements to those determined by static nucleon properties. Con- 



sider for simplicity the case of a neutral particle such as the neutron For example, the shift 
in energy due to an external electromagnetic field is determined by the zero-momentum limit 
of the amputated two-point function, 

-5M(E,B) = HmS(p) 

(T B [d-E] 2 B 2 ~E 2 2 E 2 [d 2 d-E] 

= c f9^vt~ + c D 9 + c A1 g — — c A2 g + c X3 g- 



2M ' ^ 8M 2 ' ^ 8M3 ^ 16M3 ' ^ M 4 



,2" 



B[d-E] 2 [E ■ dcr ■ B] 2 [B ■ dcr ■ E] 



c xi9 V7i 1- CxsQ r c X gg 



M 4 ; ' M 4 ; ' M 4 

2 [£V ■ dB\ 2 [B l cr ■ dE 1 } 2 cr ■ E x [d x B] . . 

+ c XW g + c xn g T7i c xl2 g — , (50) 

and we may immediately identify the various static electromagnetic moments. The identifi- 
cation of Cxi as coefficients of a well-defined effective theory is essential for correctly relating 
static coefficients, as measured by lattice [3TJ |32j |33l HH] or hadronic models, to scattering 
observables [3^1 [36J. For example, conventional definitions of the scalar electric (cxe) and 
magnetic (/?m) polarizabilities as determined by Compton scattering require that coeffi- 



cients cai and CA2 are not proportional to these quantities, but rathei^ 

4M 3 1 

a E = -c A i ~ -CA2 + Z 2 + 2Zc M + c F c s - c 2 F , 

(5 M = c A i-Z 2 . (51) 



a 
AM 3 



a 



13 Static properties of a charged particle require care in definition, see, e.g., (30] . 

14 Many approaches in the literature lack a systematic treatment of these effects [F7\ [35l [33j [3H [36] . The 



NRQED lagrangian, constrained by Lorentz invariance, trivializes relations such as ( 51 ) between static nucleon 
properties and scattering observables. 



17 



Similar relations relate static spin polarizabilities to scattering observables: see (21) and the 
associated footnote. 

The coefficients of spin-dependent operators cx% contribute to the spin-dependent structure 
functions of forward Compton scattering, obtained from the antisymmetric component of (43), 



+ H 2 (V\ iW - [f, W - [Y, 7 V) }u(k) , (52) 



where our normalization conventions are such that for a point particle MH\ = Q 2 /(Q 4 — v 2 \ 
H 2 = 0. The model- independent analysis of these functions impacts structure-dependent cor- 
rections in hydrogen spectroscopy] 1 ^] Let us compute i?i(0, Q 2 ) by using NRQED to compute 
Wjf. This requires a subclass of 1/M 5 interactions involving the magnetic field and four 
spatial derivatives, 



AC = -^-^ \ c y1 {<t ■ B, d 4 } + c Y2 {d 2 , d l a ■ Bd 1 } + c Y3 d 2 (T ■ Bd 2 + cy^t&a ■ Bd i d j 
+ c Y5 (B da 3d 2 + d 2 a- ■ dd B) + c y6 (<7 dB dd 2 + d 2 d Bcr d) 



+ c Y7 {d l B j <r ■ dd { & + &&o ■ dB J d l ) + 0(g) U 



(53) 



where we have neglected effects in setting D — > d that are not relevant to the present appli- 
cation. At low Q 2 we find 



2M Hl (0,Q 2 ) = -^Zc F + ^ 2 



Z(2c F + c s - 2c W i) - c F c D 



+ 



Q 2 



16M 4 



Cd(cwi — c f) — 2Z(c F + c W i — 32c 



Yl, 



+ 32c F (c X i + C X 2 + C X 3) + 16(-c X 7 + c X8 + cxia) 



+ 0(Q 4 ). (54) 



Using relations (15), (16), (21), and 



C Y l 



27 
256 



5 5 1 

Fx + -i-F 2 + i-F[ + ^-F' 2 + -F'{ + -F'i 



256 



16 



16 



2 ! 



(55) 



the result (54) may be expressed in terms of scattering observables. 



15 These functions enter directly in the analysis of spin-dependent transitions. Additionally, ansatzes used 
to model Wi(0, Q 2 ) may be compared to analogous studies of i?i(0, Q 2 ), which by virtue of an unsubtracted 
dispersion relation, can be reconstructed from experimental data on elastic and inelastic scattering [38 . 
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7 Summary 



We have derived a complete basis of operators and coefficient constraints through order 1/M 4 
for the parity and time-reversal invariant effective lagrangian for a heavy fermion interacting 
with an abelian gauge field, i.e., NRQED. 

The computations of this paper provide an illustration of relativity constraints in high 
orders of heavy particle effective theories jH]. In particular, the transformation law for fields 
under Lorentz boosts receives nontrivial corrections compared to a naive reparameterization 
ansatz [11]. These effects occur first in the 1/M 4 lagrangian, and are validated by the explicit 
matching and variational computations presented in this paper. Relations ^ and (25), and 



their extensions to QCD (e.g. HQET or NRQCD) represent new non-renormalization theorems 
valid to all orders in perturbation theory. 

NRQED can be used to efficiently analyze processes involving long wavelength electromag- 
netic probes of the nucleon. The analysis was motivated in part by the necessity to incorporate 
high-order corrections of proton structure in hydrogenic bound states. A sample application 
to structure-dependent two-photon exchange corrections in the small lepton mass limit was 



presented in Section 6.1 



For applications to low-energy electron scattering, m e <C E <C m p , we constructed in 



Section |5.4| the effective theory for a relativistic electron and nonrelativistic proton. The 
lagrangian and coefficient relations were derived through 0(1/M 3 ). This case may be applied 
to a model-independent analysis of radiative corrections in the extraction of proton structure 



from electron scattering, as discussed in Section |6.2| 

The effects of nucleon structure are implemented by non-pointlike values for the Wilson 
coefficients appearing in ([TJ. Having determined the structure of the NRQED lagrangian, it 
becomes a trivial task to define and compute static properties of nucleons, and unambiguously 



relate these properties to scattering measurements, as illustrated in Section |6.3 

For simplicity we focused in this paper on the effective lagrangian for a parity conserving 
theory of a heavy fermion coupled to an abelian gauge field. Extensions to nonabelian gauge 
fields, the inclusion of parity violation, and the consideration of different heavy particle spins 
each have important applications. 
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